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Abstract 

We present new micromagnetic simulations, where an influence of the 
wire shape on the spatial distribution of internal magnetoelastic stress 
is taken into account. Local magnetoelastic anisotropy is approximated 
to be uniaxial, with spatial dependence due to the stress distribution. 
The hysteresis loop is calculated with the mumax 3 software for the wire 
volume and separately for selected parts of the wire, including the surface. 
The bistability of the calculated surface hysteresis loops depends on the 
behaviour of the neighboring part of the inner core, most close to the 
surface element where the surface loop is measured. 


1 Introduction 

Amorphous micro- and nanoscopic soft magnetic wires have been found 
applications in sensors and data storage Magnetic properties 

of these wires are known to depend on the spatial distribution of internal 
and applied stresses 00. As a rule, the stresses applied are either tensile 
or torsional; these kinds of stresses come also from the glass covers and 
from the fabrication process, as it has been discussed by several authors 
0 0. On the contrary, the role of bending stresses has been largely 
ignored, with pm eed m ECU as rare exceptions. Yet, bending is known 
to strongly modify internal distribution of stresses m, and therefore the 
radius of curvature of a wire is a convenient parameter to control its 
magnetic properties. 

Our aim here is to investigate the volume and surface hysteresis loops 
of a wire in the presence of the bending stress by means of computer 
simulations. Our motivation is twofold. First is that, as suggested in m, 
the spatial distribution of stress can depend on the fabrication process. 
Therefore, stress dependent magnetic properties can vary from one series 
of samples to another. The method of simulation offers a test-bed, where 
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correlations can be investigated between theoretical stress distribution and 
the simulated magnetic properties. Second motivation is to control the 
volume and surface hysteresis loops in the same magnetization process. 
Such results are of interest when comparing and interpreting experimental 
data, obtained with the conventional induction method and with the Kerr 
effect. 

The structure of the text is as follows. In the next section we recall 
the theoretical and experimental achievements on the magnetism of bent 
wires, known in literature. In Section 3 we describe the simulation, in¬ 
cluding an approximation of an effective uniaxial anisotropy, introduced 
to meet the demands of the used software. There we describe also the 
parametrization of the spatial stress distribution, based on Refs. [10], lSj. 
Last two sections contain numerical results and their discussion, respec¬ 
tively. 


2 Literature data 

In Ref. [10], the conventional hysteresis loop was measured of a Fe-based 
wire with 125 tim diameter, as dependent on the radius of curvature R c ; 
the latter was tuned from 12 cm to about 1 cm. The bistable character 
of the loops was found to disappear for R c less than 11 cm, but it was 
restored for R c less than about 2.5 cm. The paper provides also theoretical 
evaluations of the stress distribution within the bent wire, based on the 
classical theory of elasticity m- 

In Ref. EH, the hysteresis loops have been measured by the Kerr 
method for the wires of the same diameter and composition as in Ref ESI- 
The experimental results have been accompanied by theoretical consider¬ 
ations on the energy barrier which is exceeded during the remagnetization 
process. In the discussion, the Authors highlight the role of closure do¬ 
mains and of the intrinsic stresses frozen-in during the fabrication process. 

In Ref. [12], the role of wire length is investigated experimentally for the 
same Fe-based microwires. The samples are wound on cylindrical holders, 
then the radius of curvature was equal to the radius of the holder. The 
Authors argue, that the bistability is present for any value of the radius 
of curvature, if only the wire is long enough. The study deals also with 
an influence of temperature and the amplitude of the applied field on the 
character of the loops. 

On the contrary to the literature cited above, Ref. m concerns a 
direct application of glass-covered Fe-rich microwires as bending sensors. 
The wires diameter here varies from 1 to 40 fim. Instead of the switching 
(coercive) field, the switching time is measured, what largely improves the 
sensor sensitivity. It appears that the switching time varies as the square 
root of bending. The paper contains also an updated list of applications of 
glass-coated microwires as sensors. More generally, glass-covered micro- 
and nanoscopic wires offer reduced dimensions, larger field sensitivity and 
more controllable fabrication process [5]. 
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Figure 1: The stress distribution in a straight wire. In the horizontal axis, we 
mark the position on a diameter. The wire radius r = 1 fim. 



Figure 2: The stress distribution in a straight wire, the same as in Fig. [l]- a 
cross section. White area marks the easy anisotropy axis along the wire. Colors 
online differentiate radial directions of local anisotropy axes. The positions A, 
E and F mark the wire fragments, where the hysteresis loops are calculated. 
These loops are shown in Fig. El 

3 Simulation of a Fe-rich glass-covered 
wire 

The Landau-Lifshitz equation is solved numerically by means of the mi- 
cromagnetic simulation software mumax 3 na. The diameter of the dis¬ 
cretization grains is 50 nm, while the wire diameter is 2 /xra, and the length 
is 100 iim. The Landau-Lifshitz damping constant ol— 0.1, the exchange 
stiffness constant for iron A ex = 2 x 10 -11 J/m, the magnetostriction 
constant for the Fe-rich wire A = 2.5 X 10 -6 , and the saturation magne¬ 
tization hqM — 0.6 T |5]. For the bent wire, the radius of curvature R c 
is 2 mm. The simulation is quasistatic, which is equivalent to a very low 
frequency of the field variation. 

The software mumax 3 limits the forms of magnetic anisotropy to two: 
uniaxial and cubic one. Then, numerical calculations can be performed 
only for these two forms. To fulfill this condition, we approximate the 
local magnetoelastic energy as uniaxial one. The idea is as follows. For a 
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Fe-rich magnetostrictive wire, the cylindrical symmetry allows to limit the 
number of stress components to three: radial, circumferential and axial 
one mm Locally, these three directions can be related to Cartesian 
coordinates, with three different stresses along each axis. Let us recall 
that the energy of the magnetic anisotropy is equivalent to a barrier of 
energy, which is necessary to rotate a local magnetic moment. Its uni¬ 
axial character is equivalent to the cylindrical symmetry of the angular 
dependence of the magnetic moment; if the local easy axis is along OZ di¬ 
rection, the same energy is needed to rotate the magnetic moment within 
ZX plane and ZY plane. The idea of an effective uniaxial anisotropy is 
that the magnetic moment is rotated within the plane where the energy 
cost is lower. To be more detailed, suppose that there are three tensile 
stresses a x , cr y and cr z at a given point of a magnetic sample with positive 
magnetostriction A. Without limiting the generality we can assume that 
cr x < cr y < <j z \ hence the easy magnetic axis is along OZ. The idea is 
to select the uniaxial anisotropy constant as \{a z — a y ). This choice is 
determined by the fact that the energy to rotate the magnetic moment in 
the plane ZY is less than the one to rotate it in the plane ZX. 



X 


Figure 3: The stress distribution in a bent wire. In the horizontal axis, we 
mark the position on a diameter. The wire radius r = 1 /im, and the radius of 
curvature R c =2 mm. 


In amorphous wires, the shape magnetostriction coefficient is a scalar 
m , and the spatial distribution of the magnetic anisotropy is determined 
by the distribution of stress. For the case of a straight wire, we parametrize 
the stress distribution as in nu, with small simplifications in order to 
reduce the number of different discretization grains; the latter is limited 
by the software mumax 3 to 256. The stress distribution is shown in 
Figs. Il] [2] The same information for the bent wire is shown in Figs, 
pi 4] pj These pictures are based on the calculations performed in uni, 
[ToTT In particular, the most relevant component of the bending stress a zz 
was found to vary linearly with the distance from the curvature center. 
Namely, we approximate this dependence as 


cr zz 


2 Erx 


(i) 


where E — 1.6 x 10 5 MPa is the Young modulus, r is the wire radius, R c is 
the radius of the wire curvature, and x varies within the wire volume from 
— 1 at the surface of the convex part to +1 at the surface at the concave 
part. The approximation in Eq. [l] is justified as long as the wire radius r 
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Figure 4: The stress distribution in a bent wire - a cross section. The distri¬ 
bution is as in Fig. [3j with the orientation of the stress along OX (for details 
see the text). White area marks the easy anisotropy axis along the wire. The 
positions marked with labels indicate the wire fragments, where the hysteresis 
loops are calculated, as shown in Figs. [l2]and[T3j The color online on the right 
part of the figure is related to the assumed direction (but not the value) of the 
magnetic anisotropy in the plane perpendicular to the wire. 


is much less than the radius of curvature R c . For r/R c — 5 X 10 4 , the 
bending component of the stress changes from -800 to +800 MPa. 

Basically, we expect that in the compressed part of the wire a competi¬ 
tion exists between the stresses aee and a rr ■ This competition leads to a 
different direction of effective stress at each point of the wire. This fact is 
unfavourable, as we have to keep in computer memory each grain of dis¬ 
cretization. To evade this difficulty, we assume a homogeneous direction 
of the magnetic anisotropy (but not the value) in the compressed part of 
the inner core. These two versions of the model distribution of stress are 
shown in Figs. am To check this assumption, calculations on the bent 
wire are performed twice: for this direction either OX (Fig. [4]) or OY 
(Fig. Below we show that the results do not differ remarkably; hence 
the assumption is justified. 


4 Results 

In Fig. [6] we show the spatial distribution of the magnetization in a 
cross section of a straight wire of radius r = 1/xra, at half of its length 
L [L — 100 iim). Both the outer shell, with domains oriented radially, 
and the inner core - one large domain, uniformly magnetized along the 
wire, are visible. This result reproduces what is known about the domain 
structure of wires with positive magnetostriction A [HI [20j 112- 

In Fig. [7| the distribution of magnetization is shown for similar cross 
sections at small but different distances from the wire end. From these 
pictures, it is clear that the axial component of the magnetization increases 
with the distance from the wire end. This is consistent with what is known 
on the closure domains in Fe-rich wires [221120] . Also, we have calculated 
the hysteresis loops for wires of different lengths. These results are shown 
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Figure 5: The stress distribution in a bent wire - a cross section. The distribu¬ 
tion is as in Fig. [3j with the orientation of the stress along OY (for details see 
the text). White area marks the easy anisotropy axis along the wire. The color 
online on the right part of the figure is related to the assumed direction (but not 
the value) of the magnetic anisotropy in the plane perpendicular to the wire. 



Figure 6: The magnetization distribution in a cross section of a straight wire, at 
the half of the wire length. The colors (red and blue online) mark the horizontal 
component of the orientation to right and left, respectively. 


in Fig. [8j and the remanence as dependent on the wire length - in Fig. 
[9] These results show that the bistability emerges for the wire length 
larger than some critical value, as in experiment ESI- Yet the values of 
this critical value depends on the wire diameter, and there is no one-to- 
one correspondence in the results of the calculations and the experimental 
data. 

A test of the approximation of two anisotropy orientations in the com¬ 
pressed volume is to compare the hysteresis loops for both these cases. 
We show the results of this test in Fig. |10| as a rule, only the upper 
branch of the loop is shown. In our opinion, the results are satisfactory; 
two curves approximately coincide and none of them is bistable. On the 
contrary, the bistability is present for the straight wire. 
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Figure 7: The magnetization distribution in a cross section of a straight wire at 
the distances zero, 0.25 and 0.5 /im from the wire end from top to bottom. The 
colors (red and blue online) mark the horizontal component of the orientation 
to right and left, respectively. 
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In Fig. [lT]we show the volume and surface hysteresis loops for a straight 
wire. Basically, in the lack of bending stress the three surface loops, taken 
from three sides, should be the same. Yet, some quantitative differences 
are visible, which indicate an inaccuracy of the numerical method. 

In Figs. |12|and p~3] we show the hysteresis loops which characterize the 
magnetic state of some fragments of the wire, on the surface and within 
the wire volume. The surface loops can be measured by the Kerr method. 
The internal fragments cannot be measured without removing the surface; 
however, such removal modifies the stress distribution. Therefore the 
simulation is the only way of an insight into the magnetic behavior of 
internal parts of the wire. 


The results shown in Fig. [12] indicate, that the fragments in the left 
part of the wire, where the bending stress is tensile, produce bistable 
contributions to the hysteresis loop. For first two fragments from the left 
(in the positions A and B in Fig. [4|, the loops almost coincide, with 
the coercive field about 11 kA/m. Just on the other side of the wire 
center, where the bending stress becomes compressive, the coercive field 
is abruptly reduced (position C in Fig. [4] dark blue line online in Fig. 


12). Quite surprisingly, a jump of magnetization is seen as far as in point 


D; it is only the remanence what is reduced (red line online in Fig. 12) 


On the contrary, at the surface (position E in Fig. [4] blue line online i 
Fig. 12)) there is no hysteresis at all; this is because at this area, the local 
magnetic anisotropy is radial and strong. 



B z [kA/m] 


Figure 8: The upper branches of the hysteresis loops for a straight wire, for 
different wire lengths. 


In Fig. |13| the contributions to the hysteresis loop are shown for the 
fragments F, G, and H in the Fig. [4] At the position F, the bending 
stress is neither tensile, nor compressive; its planar contribution is ne¬ 
glected here. The related hysteresis (red line online in Fig. 13) is far from 


bistability. On the contrary, the contribution at point G (within the inner 
core) shows a ’partial bistability’; apparently, the fragment G consists of 
two antiparallel domains. Here again we see a difference between the left 
and right part of the wire. A similar behavior is seen at position H, right 
in the wire center (dark blue online in Fig. 13), with some departure from 
the preceding case near the coercive field. 
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Figure 9: The remanence magnetization for a straight wire, as dependent on 
the wire length - the latter in logarithmic scale. The transition is visible from 
non-bistable to bistable behavior, when the wire length exceeds some critical 
value. 



Figure 10: The upper branch of the hysteresis loop in a straight wire (red online) 
and of a bent wire, with two orientations of the easy axis in the compressed part 
of the wire: OX and OY (green and blue online). The wire radius r = 1 /im, and 
the length L = 100 /xm. For the bent wire, the radius of curvature R c = 2 mm. 

5 Discussion 

In our simulations, two assumptions have been induced by the demands of 
the software mumax 3 and the aim to take into account details of the do¬ 
main structure in a bent wire. First was the approximation of an effective 
uniaxial anisotropy, and second - a homogeneous direction of the mag¬ 
netic anisotropy perpendicular to the wire axis in its compressed part. 

The quality of the first assumption has been verified by the proper re¬ 
production of the domain structure in the straight wire. To be more 
precise, this assumption neglects the difference between the actual angu¬ 
lar dependence of the energy of the magnetization and the model uniaxial 
anisotropy where we have just an easy axis. In the latter model, a local 
magnetic moment can deviate from its easy axis in any direction with the 
same energy cost, while in real situation a rotation within one plane is 
easier than within another plane. Yet, the approximation of an effective 
uniaxial anisotropy does not destroy the domain structure, as seen in Fig. 

m 
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Figure 11: The hysteresis loops, volume (red online) and three surface ones 
(blue, green and violet online) in a straight wire. The wire radius r = 1/ira, the 
length L = 100, and the radius of curvature R c = oo. 


Our second assumption is that in the compressed part of the wire, the 
actual direction of the magnetic anisotropy perpendicular to the wire is 
not relevant. If this is true, we can set this direction along an arbitrary 
axis without modification of the final results. Actually we have performed 
this test, as shown in Fig. [To] There, the two curves for the bent wire 
almost coincide. This means that the assumption does not influence the 
results qualitatively. Also, our approach is validated by reproduction of 
two experimental effects: the closure domains and the critical length for 
the bistability. 



B z [kA/m] 


Figure 12: The hysteresis loops - contributions from different parts of the wire, 
from a piece of surface at the concave side (bistable loops, red and green online), 
to just after the center (dark blue online), to near the surface on the convex, 
compressed side (violet online) and at the surface (blue online, no hysteresis). 
The wire radius r = l//m, the length L = 100, and the radius of curvature 
R c = 2 mm. 


The results of our calculations performed for the bent wire indicate, 
that below some value of R c , the bistability is destroyed. Here we are not 
in position to state if the bistability is restored if the radius of curvature 
is even less, as it was found experimentally in M- New result of our 
calculations is that the stress inhomogeneities inside the bent wire produce 
a specific domain structure. As for straight wires, this structure is the 
result of an interplay between the magnetoelastic anisotropy, the exchange 


10 











0.5 



-12 -10 -8 -6 -4 -2 0 2 4 

B z [kA/m] 


Figure 13: The hysteresis loops - contributions from different parts of the wire, 
from a piece of surface F as shown in Fig. [4] (red line online), up to position G 
(green line online) to the center, position H (dark blue online). The wire radius 
r = 1 /xm, the length L = 100, and the radius of curvature R c = 2 mm. 


and the pseudodipolar interaction. However, in a bent wire the cylindrical 
symmetry is broken by the bending stress. As a consequence, the inner 
core is shifted to the convex part, and the compressive bending stress in 
the concave part of the wire produces a strong anisotropy perpendicular 
to the wire axis, which destroys the bistability. We conclude that bent 
wires offer more options of spatial configurations, as ordered and random 
arrays [241125] . and therefore enable new opportunities of applications. 
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